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ABSTRACT 



Star formation in galaxies is observed to be associated with gamma-ray emis- 
sion, presumably from non-thermal processes connected to the acceleration of 
cosmic-ray nuclei and electrons. The detection of gamma rays from starburst 
galaxies by the Fermi Large Area Telescope (LAT) has allowed the determination 
of a functional relationship between star formation rate and gamma-ray luminos- 
ity ( lAckermann et al.l 120121 ) . Since star formation is known to scale with total 
infrared (8 — 1000 jum) and radio (1.4 GHz) luminosity, the observed infrared and 
radio emission from a star-forming galaxy can be used to quantitatively infer the 
galaxy's gamma-ray luminosity. Similarly, star forming galaxies within galaxy 
clusters allow us to derive lower limits on the gamma-ray emission from clusters, 
which have not yet been conclusively detected in gamma rays. In this study we 
apply the functional relati onships between g a mma -ray luminosity and radio and 



IR luminosities derived in Ackermann et al 



(120121 ) to a sample of the be s t can- 



didate galaxy clusters for detection in gamma rays from lAckermann et al.l (120101 ) 
in order to place lower limits on the gamma-ray emission associated with star 
formation alone in galaxy clusters. We find that several clusters have predicted 
gamma-ray emission from st ar formation that are wi thin an order of magnitude 
of the upper limits derived in lAckermann et al.l (120101 ) based on non-detection by 
Fermi-LAT. Given the current gamma-ray limits, star formation likely plays a 
significant role in the gamma-ray emission in some clusters, especially those with 
cool cores. We predict that both Fermi-LAT over the course of its lifetime and 
the future Cherenkov Telescope Array will be able to detect gamma-ray emission 
from star-forming galaxies in clusters. 
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INTRODUCTION 



Galaxy clusters are the most massive gravitationally bound objects in the universe, 
and contain a dense population of galaxies surrounded by gas distributed throughout the 
intracluster medium (ICM). The complex environment of galaxy clusters is bound to host, 
at some luminosity level, processes that lead to the production of gamma rays. In particular, 
a significant fraction of the gamma-ray emission is thought to be associated with the ICM, 
and more specifically with cosmic ray (CR) populations accelerated by shocks and turbulence 
wi thin the ICM , as well as, possibly, with dark matter annihilation and decay (first suggested 



m 



Tot anil |2 004 ) . In the ICM, CR protons and other nuclei collide with dust and gas particles 



and decay to charged and neutral pions, which in turn decay into leptons and gamma rays 
respectively. CR electrons, if sufficiently energetic, can inverse- Compton s catter photons 



from the Cosmic Microwave Background up to gamma-ray energies (see e.g.. iRephaeli et al. 



20081 for a general description of nonthermal activity in the ICM of clusters of galaxies). 



Clusters have not, however, been detected yet as gamma-ray sources. The Large Area 
Telescope (LAT), the primary instrument aboard the Fermi Gamma-ray Space Telescope, 
has placed upper limits on the flux from the best candidate cluste rs and used these uppe r 
limits to place constraints on the cosmic ray populations in clusters (jAckermann et al.ll2010l ). 
Atmospheric Cherenkov telescopes such as H.E.S.S. and MAGIC hav e also reported nu ll 
results from observations of sele cted clusters of galaxies, s uch as Perseus (lAleksic et al.ll2012l) . 
Coma flAharonian et al.ll2009ah . and A0085 and A0496 (Uharonian et al.ll2009bh . 



In addition to the gamma rays produced by the ICM, the cluster member-galaxies 
themselves are at some level a source of gamma rays. Ordinary galaxies such as our own 
Milky Way, its satellite galaxies the Large Magellanic Cloud (LMC) and Small Magel- 
lanic Cloud (SMC), and the nearby Andro meda Galaxy (M31) have a ll been detected in 



gamm a rays (MW: lAbdo et all |2009j , LMC: lAbdo et all l2010dl SMC: lAbdo et al.l 12010b 
M31: lAbdo et al.l l2010d ). In additio n, the recent detec tion by the Fermi-LAT of four 



star-forming galaxi es, M82, NGC 253 (lAbdo et al.l l2010ah and NGC 1068 and NGC 4945 
( jLenain et al.ll2010l ). allowed for the determination of a quantita tive functional relation ship 
between the star formation rate and the gamma-ray luminosity ( jAckermann et al.l 120121 ) . 



While gamma-ray emission has not yet been detected from galaxy clusters, some min- 
imum gamma-ray emission resulting from star formation activity in cluster members must 
exist. Lower limits on the gamma-ray flux from clusters of galaxies can therefore be deter- 
mined by considering only the emission from cluster members with ongoing star formation, 
and can also provide insight into the star-forming population of galaxies within clusters, 
as compared to the field. This is the scope of the present study. The relations derived 



in 



Ackermann et al.l ( 120121 ) are used to predict lower limits on the gamma-ray luminosity 
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coming from star formation within cluster members alone, provided total IR and/or radio 
luminosity measurements for a sample of clusters. We then compare thes e estimates to the 



Ferm i upper limits on the gamma-ray luminosity from the same clusters (lAckermann et al. 



20101 ) . and assess the potential of the Fermi-LAT over its lifetime as well as of the future 
Cherenkov Telescope Array for detection of gamma-rays from star-forming galaxies in clus- 
ters. 

In the following section, we review the sources of multiwavelength emission resulting 
from star formation in galaxies and clusters. In section 3, we describe our cluster sample 
selection and the available IR and radio data. In section 4 we describe our results and present 
lower limits on the gamma-ray emission from a selection of clusters. We use the lower limits 
to explore the possibility of detection by various gamma-ray telescopes in section 5, and 
conclude in section 6. 



2. EMISSION FROM STAR FORMATION 

The bulk of the bolometric luminosity of young stars is emitted in the ultraviolet (UV), 
which is then efficiently reprocessed into infrared (IR) light by dust in the interstellar medium 
(ISM) within galaxies su ch that their spectral energy distributions are peaked in the IR 



(jSanders fc Mirabel! Il996l ). The total IR luminosity (8 — 1000/xm) thus traces closely the 



rate of star formation in a galaxy and is an extinction-fr ee way to measu re the star for- 



mation rate (SFR) since the ISM is optically thin to IR (IKennicuttl Il998l ). However, the 



total IR luminosity is not directly measured; rather, single-band measurements (24 /im and 



Dale & Helou 


20021) or using r 


1 (e.g., 


Sanders & Mirabel 


1996) 



modeling the dust emission in the ISM 
common in the literature and a determination of the total IR luminosity thus depends some- 
what on the particular templates or relations u s ed fo r con version. Discrepanc i es be tween 
methods developed specifically by lDale fc Helou! (120021 ) and lSanders fc Mirabel! (Il996[ ) vary 
from ~ 4% for starburst galaxies, which have high SFRs and ge nerally are very IR - bright 
(L 8 _ioooMm > lO n L ), to ~ 25% for low-luminosity galaxies, with lSanders fc Mirabe 
method predicting smaller values for L 8 _i oo/mi 



than those predicted by 



Dale fc Helou 



(1996) 

SI). 



Massive stars typically end their life-cycles as sup ernovae, whose remnants (SNRs) are 



likely accelerators of cosmic ray protons and electrons (e.g. , IVolklll989l and references therein) . 
The CR electrons interact with magnetic fields within a galaxy and quickly lose their energy 
through synchrotron radiation in the radio continuum. Measured at 1.4 GHz, the radio 
luminosity is tightly correlated with total IR luminosity in star-forming galaxies over sev- 
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eral orders of m agnitude and is therefore also used as an indicator for SFR ( jCondonlll992 
Yun et al.ll200ll ). This correlation can be explained by a si mple mode l that treats galaxies 
as "calorimeters" of both UV photons and CR electrons (jVork! Il989l ). CR electrons lose 
most of their energy through synchrotron radiation and UV photons lose most of their en- 
ergy through absorption and re-emission by dust in the ISM. In this model, if the ratio of 
synchrotron emission versus total energy losses in the CR electron population is large, then 
the radio luminosity becomes a proxy for SFR. In addition, if UV photons emitted by stars 
are effectively reprocessed into the IR, then IR luminosity would also trace with SFR. The 
Li.4GHz-L 8 -ioooMm ratio, which should be constant if both are due mainly to star formation, 
can be thus be used as an independent indicator of the SFR in galaxies. Active galactic 
nuclei (AGN), which can have large luminosities in both IR and radio due to processes other 
than just star formation, and which presumably also con tribute to the ga mma-ray emission 
in some clusters, deviate from this IR-radio correlation (lYun et al.ll200ll ). We discuss this 
caveat in detail for the cluster sample under consideration here. 

Gamma rays are expected to scale with SFR in normal, star-forming galaxies, and are 
thought to arise mostly from CR nucleon collisions with the ISM and subsequent decay into 
charged and neutral pions, which then decay to leptons and gamma rays, respectively. At 
energy thresholds higher than those required for pion decay the gamma-ray spectrum is 
the same shape as the underlying CR nuclei population, which is controlled by the number 
of SNRs. Therefore the amplitude of the gamma-ray spectrum can be used as a proxy for 
SFR. This is supported by observations of the LMC in gamma rays, which has been spatially 
resolved by Fermi-LAT. The LMC's star forming region 30 Doradus is associated with strong 
diffuse gamma-ray emissio n, which implies the CR proton intensity is also strongest in that 
region (jAbdo et al.ll2010dl ). The secondary CR electron population, resulting from charged 
pion decay, also contri butes to the gamma-ray spectrum primarily by inverse- Compton up- 
scattering of starlight (lAckermann et al.ll2012l ). 



Early in the miss ion, in all-sky surve y mode, Fermi-LAT detected two starburst galaxies, 
M82 and NGC 253 (jAbdo et al.l l2010al ). Those starbursts were also detecte d in the very- 



high- energy regime by ground-b ased gamma-ray te lescopes (M82 by VE RITAS; I Acciari et al. 
20091 and N GC 253 by H.E.S.S . : lAcero et al.l l2009h. As first reported in lLenain et all (bold ) 



and later in lAckermann et al.l (120121 ). F ermi-LAT detected two a dditional starbursts, NGC 
1068 and NGC 4945, which host AGN. lAckermann et al.l (120121 ) examined a sample of 69 
galaxies, mostly consisting of starbursts for which Fermi-LAT has upper limits, 9 of which 
are associated with AGN. The sample includes the four starbursts detected by Fermi-LAT 
and several detected local group galaxies which are not starbursts: th e Sma ll and Large 
Magellanic Clouds, the Milky Way, M31 and M33. lAckermann et all (boij ) found that 



gamma-ray luminosity scales with total IR and radio luminosities and hence with SFR across 
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several orders of magnitude in luminosity and SFR. The relationships found between SFR 
and gamma-ray luminosity are power-law in nature and differ slightly depending on whether 
AGN are included in the sample, but the dispersion in the relationships dominate over the 
differences in the fits. 

Star formation in cluster galaxies occurs at slower rates compared to that in field galaxies 



^see e.g.. lReddy fc Yunll2004l . and references therein ]. However, the pro cesses that suppress 
star formation in clusters are not well-understood. IChung et al.l (120111 ) used infrared data 
from WISE (Wide-field Infrared Survey Explorer) to estimate SFRs for a sample of local 
clusters with known masses and found no correlation between cluster mass and SFR. This 
implies that clusters processes that scale with mass perhaps only weakly affect star formation, 
if at all. Ram pressure in particular scales up with cluster mass and is thought to influence 
star formation. Ram pressure strips gas from a galaxy as it moves through the ICM, removing 
the raw material needed for star formation. Rates of star format ion in cluster gala xies 
get smaller as the radial distance from the cluster center decreases (j Chung et al.l 1201 ll and 
references therein) which supports the theory that ram pressure stripping strongly influences 
star formati on. In addition, ram p ressure is shown to decrease star formatio n in clusters in 



simul ations (IBook fc Benson! 120101 ) and in observations of the Virgo cluster (jVollmer et al. 



20081 ). However, this conflicts with the lack of evidence for a corre lation between mas s and 



star formation in local clusters, described by IChung et al.l ( 1201 ll ) and others (e.g., iGoto 



20051 ). Galaxy-galaxy mergers and interactions also tend to trigger bur sts of star formation 



in clu sters, but are most common in poor clusters and groups, (e.g., iMartig &: Bournaud 
20081 ). while the best clusters for detection in gamma rays tend to be rich and massive. 
Still many cluster galaxies do show evidence of ongoing star formation, and the density of 
galaxies within clusters mean that the summed contribution to the gamma-ray emission may 
be significant. 

Cool-core clu sters display elevated S FRs in their central galaxies, as compared to non- 
cool-core clusters ( IMcDonald et al.ll201ll ) . These clusters have high density, low temp erature 
cores which suggest that cool gas is steadily flowing into the center (for a review, see iFabian 
19941 ). Therefore the elevated SFR may be due to this mass deposition onto the central 
galaxy, which would enhance star forma tion by supply i ng th e brightest cluster galaxy (BCG) 
with a steady stream of raw material. lEgami et al.l (120061 ). however, found only a weakly 
suggestive relationship between mass deposition rate and SFR in a sample of X-ray selected 
clusters. In addition to several clusters of galaxies, in this work we also consider the possible 
gamma-ray emission from star formation in 3 BCGs located in cool-core clusters. 
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METHODS 



We estimate here the gamma-ray luminosity associated with star formation in cluster 
galaxies from the total infrared luminosity (8 — 1000 /xm) and from the 1.4 GHz radio contin- 
uum luminosity using the relationships in lAckermann et al.l ( 120121 ). We restrict our analysis 
to emission from individual galaxies where there is active star formation, since star forma- 
tion in clusters only occurs within galaxies. We performed a lit erature search for I R an d 
radio data for cluster galaxies starting with the list of clusters in lAckermann et al.l ( 120101 ). 
which places upper limits on the gamma-ra y flux for the best cand idate clusters based on 
nondetection by Fermi-L AT . The clusters in Ackermann et al. ( 20icf) were se lected from the 
HIFLUGCS catalog of brightest X-ray clusters ( jReiprich fc Bohringerl 120021 ). The clusters 
with the highest mass-to-distance-squared ratios were selected for analy sis. Additionally, 



several clusters with nonthermal radio emission were included in the list of Ackermann et al. 



( 120101 ). We performed a literature search for radio and IR data on all of the clusters in that 
list, excluding several that were close to the galactic plane or that had high redshifts (with 
the exception of the Bullet Cluster). We excluded AGN from our calculated luminosities, 
since AGN can be bright in IR and radio but the emission is not due to active star forma- 
tion. The results of the literature search are presented in Table 1 and the corresponding 
estimated gamma - ray lu minosities are summarized in Table 2; the upper limits provided by 
Ackermann et al.l ( 120101 ) are also included in Table 2. 



In looking for IR and radio data, we preferentially searched first for total IR (8 — 1000 
jum) and 1.4 GHz luminosity function s (LFs) of cluste rs with Schechter fits. W e found fitted 
total IR LFs for two clusters, Coma dBai et al.l 120061) and the Bullet Cluster (jChung et al. 
2010h . While we fou nd a radio LF for Coma jMiller et aD l2009h and several far-IR (100-500 
jum) LFs for Virgo (IDavies et al.l 120101 ). they were not fitted and therefore we were unable 
to convert to total luminosities. 

We expanded th e search to include lists of luminosities of individual clust er members. 



Reddy fc Yunl (120041 ) presented L L4GHz and L 60/im for members in 7 clusters. iRieke et al. 



(120091 ) found a tight correlation between L 60Aim and Lg^xooo^m and derived a r elationship for 
conver ting one to the other; we used this relation to convert the 60 /im data in lReddy &: Yun 
(120041 ) to total I R. We converted each rad io and IR galaxy luminosity to gamma-ray using the 
relationships in lAckermann et al.l ( 120121 ) that exclude AGN, then summed up the gamma- 
ray luminosities of individual galaxies to obtain a total cluster gamma-ray luminosity. We 
excluded galaxies classified as Seyferts or LINERs, which are typically AGN, from our cal- 
culations. Any galaxies not identified as Seyferts or LINERs that had one anomalously large 
luminosity (radio or IR) as compared to the other were checked in the l iterature, and those 
identified as hosting an AGN were also excluded. In lReddy fc Yunl (120041 ). only cluster mem- 



-7- 



bers detected in both radio and IR were included, so total cluster luminosities using that 
data are most likely underestimates in both IR and radio. In the cases for which we had lu- 
minosity functions, we first integrated the LF, then converted this total cluster luminosity to 
gamma-ray. Since the relationship between IR and gamma-ray luminosity is slightly steeper 
than linear, this procedure results in an overestimate of the true total gamma-ray luminosity 
by ~ 20%. This, however, is smaller than the intrinsic dispersion in the relationship reported 



m lAckermann et al.l ( 120121 ). which we take as the uncertainty in our estimates of gamma-ray 
luminosities. 

Uncertainties in galaxy IR and radio luminosities were not reported in the papers we 
found. We therefore report the uncertainties in the estimated gamma - ray lu minosities as the 
dispersion in the power-law relationships derived in lAckermann et al.l (120121 ) . For L .i_iooGeV 
estimated from L s _ 1000fim , the uncertainty is 0.25 dex (excluding AGN ). For L n .i-innGev 
estima ted from Li.4ghz, the uncertainty is 0.19 dex (excluding AGN). While lAckermann et al. 
(120121 ) provides uncertainties in the fit parameters, the scatter in the Lo.i-iooGeV-Ls-iooo/xm 
and Lo.i_iooGeV-Li.4GHz relationships dominate over the fit parameter uncertainties. 



In coo l core clusters, there is evidence for active star formation in the central cluster 



galaxies ( IMcDonald et al.l 1201 ll ). We found radio fluxes for the brightest cluster galaxy that 
are not associated with AGN in 3 cool core clusters, Ophiuchus, A2029, and A2142, and 
included them in our sample. While Ophiuchus is close to the galactic plane, we include 
it as a representative cool-core cluster. The gamma-ray luminosities estimated from these 
clusters are most likely severe underestimates, but it is interesting to note that for cool core 
clusters the BCG alone can have a significant predicted gamma-ray emission due to increased 
star formation. 

We note that there is observed large scale (~ 1 Mpc) diffuse radio emission from the 
ICM of some galaxy clusters which indicate the existence of cosmic rays tha t may also be 
responsible for as-yet-undetected gamma-ray emission. (IFerrari et al.l 120081 ) . These ICM 
cosmic rays are accelerated primarily by shocks and turbulence in the ICM, as opposed to 
supernova remnants in cluster galaxies, and are not related to star formation in galaxies. 
Observations of clu sters in the IR are always associated with individual galaxies (see e.g., 
Coppin et al.ll201ll ) . We thus choose to consider only the radio and IR emission from galaxies 
themselves in order to place conservative lower limits on the gamma-ray emission from 
clusters from star formation in cluster galaxies alone. 
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RESULTS 



We present our results in Table 1 and 2 and in the associated Figures 1 and 2. Table 
1 shows the IR and radio luminosities for the clusters under consideration, while Table 2 
presents our calculation for the gamma-ray luminos ity lower limits, and compares the latter 
with the upper limits from lAckermann et al.l ( 120101 ) . Finally, Figures 1 and 2 show the cal- 
culated gamma-ray luminosities as a function of the IR and radio luminosities, respectively. 

For several clusters the lower limits on gamma-ray luminosity predicte d from star forma- 



tion a re within an order of magnitude or so of the upper limits derived in lAckermann et al. 



(120101 ) . The lower limit luminosity of the BCG (IC 1101) in the cool-core cluster A2029 in 



imit. The core of A2029 shows extended radio emis- 



particular is quite close to the upper 
sion which suggests it hosts an AGN (ITaylor et al.lll994l ). A2029 has also been observed in 



X-rays, but it is unclear whether the X-ray emission is due to an AGN since the emission is 
not point-source-lik e and does not fit a power-law spectrum, as is typical for X-ray AGNs 
(IClarke et al.l 120041 ) . Additio nally, the star formation rate calculated nominall y from the 



presen ted radio luminosity in iMcDonald et al.l ( 120111 ) using the relationship in iKennicutt 



(119981 ) yields an SFR that is higher than that of any local starburst galaxy; therefore is it 
most probable that only a fraction of the radio luminosity is due to star formation. However, 
the large ratio of far UV to Ha emission sugge sts that the BCG of A2 029 is an older starburst 
galaxy with ongoing star formation activity ( IMcDonald et al.l 1201 if ), so it is included in our 
sample. 

The Bullet Cluster may not be detectable by Fermi-L AT due to its distance (z ~ 0.3). 
It is however an interesting cluster, with a total IR LF available in the literature and a high 
IR luminosity, and was therefore included in our sample. If a cluster similar in mass and 
activity to the Bullet existed nearby it would likely be observed in gamma rays. The Bullet 
Cluster's predicted emission from star formation is much lower than the upper limit on its 
gamma-ray luminosity given its distance, but this lower limit on the Bullet's gamma-ray 
luminosity is comparable to the detection threshold for a closer cluster with similar mass 



such as Coma (although the IR LF is a different shape than that of Coma ( iChung et al. 
2O10h ). 



We present two tota 
IR luminosit y function ( 



IR luminosities for the Coma cluster, one calculated from the total 



Bai et al.l 120061 ) and the other from a sum of individua l member 



l umin osities ( IReddy &: Yunll2004l ). The galaxy sample used to calculate the LF in lBai et al. 



( 120061 ) does include 3 AGN, but their contribution to th e total luminosity is negligible. 
We therefore used the relation in lAckermann et al.l ( 120121 ) that excludes AGN, as for all 
other clusters. As expected, the gamma-ray luminosity as calculated from the IR luminosity 
function is larger than the gamma-ray luminosity calculated from summing individual galaxy 
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luminosities. 



DISCUSSION 



5.1. Detections with Fermi-LAT 



Given the lower limits presented in the previous section, we examine here which clusters, 
if any, may be detectable by Fermi-LAT over the course of the instrument's lifetime. Fermi- 
LAT is an all-sky survey telescope, covering the full sky every 2 orbits, or 3 hours. It 
was launched in June 2008 and is currently funded through at least 2016; it is expected 
to be operational for a total of 10 years. Its sensitivity gets better over time, and if we 
assume a 10-year flux sensitivity of ~ 3 x 10~ 9 cm -2 s -1 for detecting a point source with 
phot on index = 2.2 at the 5a level, which is the best fit value for star-forming galaxies 
from lAckermann et al.l (120121 ) , then several clusters have lower limit fluxes that are above 
or within a factor of 2-3 of this sensitivity limit, including A2029, Virgo, Coma, A1367, and 
Hydra, with the flux of A2029 strictly over the 3 x 10~ 9 cm 



'' 2 s 1 l imit. The photon index 



a 1 = 2.2 is chosen as it is the value used in lAckermann et al.l ( 120121 ) to predict upper limits 
on gamma-ray luminosities from non-detected starburst galaxies. However, Virgo is a large, 
extended object, and it hosts a bright AGN, M87, which ma y make it d ifficul t to detect 
emission in excess of this gamma-ray point source. Recently, lHan et al.l ( 120121 ) claimed a 
detection of diffuse gamma-ray emission from the central 3° region of the Virgo Cluster due 
to dark matter annihilation; however further analysis is needed to confirm or deny this claim. 
Additionally, the BCG in A2029 may also host an AGN, which would also make detection of 
gamma-ray emission from the cluster as a whole potentially difficult. Coma specifically is a 
good candidate for detection. Rec ent radio observations of Coma show diffuse radio emission 
in the form of a halo and a relic (IBrown fc Rudnickll201ll ). indicating the existence of non- 
thermal processes throughout the cluster, which imply the existence of gamma rays, and as 
shown here star formation in the galaxies within Coma can also give significant gamma-ray 
emission. 



5.2. Prospects for Detection with Ground-based Gamma-ray Telescopes 



Galaxy clusters have not yet been detected by ground based gamma-ray telescopes, 
such as H.E.S.S., MAGIC or VERITAS. Recent studies of the Perseus cluster by MAGIC 
did not yield a detectio n apart from the central AGN, NGC 1275 and a radio galaxy, IC 
310 (lAleksic et al .1120121 ). However, the gamma-ray spectrum of NGC 1275 as measured by 
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MAGIC drops off above 630 GeV, which means that Perseus is a good potential candidate 
for detection in this energy range. The lower limit gamma-ray flux of Perseus above a 
TeV, about 1 x 10~ 19 c m" 2 s -1 , is approxi mately 6 orders of magnitude smaller than the 
upper limits reported in lAleksic et al.l ( 120121 ). and well below the sensitivity limits of current 
atmospheric cherenkov telescopes. The gamma-ray emission associated with star formation 
in Perseus is therefore unlikely to be detected by ground-based instruments. 

We also investigated whether the next-generation ground-based telescopes, such as the 
Cherenkov Tele scope Array (CTA ), would be able to produce a detection of galaxy clusters. 
As described in lActis et al.l (j201l[ ). the effective area of CTA is typically the limiting factor 
for a given observing time, usually about 25-50 hours. If we assume a constant power law 
spectrum for the gamma-ray emission from star-forming galaxies in clusters with a photon 
index a 7 = 2.2, we can estimate the lower limits on the flux from the clusters in our sample 
from 100 GeV to 10 TeV, which is the primary target energy range for CTA. The cluster 
with the largest flux in this range is A2029 with 2.6 x 10 -15 cm~ 2 s -1 ; Virgo then Coma have 



the next highest fluxes at 6.7 x 10 



-16 cm 2 g i anc j x ]_q 16 cm 



-2 



s 1 , respectively. A1367 



and Hydra also have predicted fluxes that are a factor of ~ 2 lower than Coma's. Using the 
flux of Virgo, the required effective area of CTA for 50-hour observation would be about 8 
km 2 ; CTA will ideally cover at tens of km 2 , at least. These est imates assume no spectral 
breaks at higher energies; as described in lAckermann et al.l (120121 ) . starburst galaxies appear 
to have a spectrum described by a single power-law but the spectra of local galaxies such 
as the Milky Way are better described with an exponential cutoff or broken power-law. We 
therefore take these very-high-energy fluxes to be optimistic predictions. We conclude that 
the performance of CTA may allow for the detection of gamma rays from star formation in 
galaxy clusters for potentially several clusters. 



CONCLUSIONS 



We calculated lower limits on the gamma-ray emission from galaxy clusters consider- 
ing only cluster member galaxies with active star formation using observed IR and radio 
luminosities for s e lected , nearby massive clusters. Employing the relationships derived in 
Ackermann et al.l ( 120121 ) for Fermi-LAT detected galaxies, we converted IR and radio clus- 
ter luminosities into gamma-ray luminosities. Several clusters have lower limits on their 
gamma-ray emission that are within ab out an order of magnitud e of the upper limits based 
on the Fermi-LAT non-detections from lAckermann et al.l ( 120101 ). implying that star forma- 
tion could contribute at the level of 10% or more to cluster gamma-ray emission. Several 
clusters also have lower limits that are within a factor of a few of the 10-year sensitivity 
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limits of the Fermi-LAT; the best candidate clusters for detection by Fermi-LAT based on 
these lower limits are A2029, Virgo, Coma, A1367, and Hydra. CTA will also likely be able 
to detect clusters such as A2029, Virgo and Coma for anticipated instrumental performance 
and design. Star formation may thus be a significant source of gamma-ray emission for some 
galaxy clusters. 

This work is partly supported by NASA grants NNX09AT96G and NNX09AT83G. SP 
acknowledges support from an Outstanding Junior Investigator Award from the Department 
of Energy, DE-FG02-04ER41286. 
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Table 1 
IR and Radio Luminosities 



Name 


l°g(L8-1000Atm) 

(L ) 


Ref. 


log(L L4G Hz) 

(W Hz" 1 ) 


Ref. 


D L 

(Mpc) 


Coma (LF) a 


11.75 


2 






100 


Coma 


11.32 


1 


23.16 


1 


100 


AWM7 


10.80 


1 


22.28 


1 


69.2 


Perseus 


10.53 


1 


22.05 


1 


72.3 


Hydra 


11.10 


1 


22.68 


1 


57.1 


A1367 


11.48 


1 


23.25 


1 


96.6 


Virgo 


10.86 


1 


22.42 


1 


19.4 


Bullet 


12.63 


3 






1479 


Ophiuchus 6 






22.68 


4 


118 


A2029 6 






24.86 


4 


339 


A2142 6 ' c 






22.68 


4 


401 



Notes. Dl is luminosity distance, retrieved from the NASA/IPAC Extragalactic Database 
(NED). The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 

( a ) LF is luminosity function. See text for details. 

( b ) The data for these clusters is for BCGs only, not the full clusters. 

( c ) The radio lum inosity for the BCG o f A21 42 is an upper li mit- 
References. (1 ) iReddv fc Yunl J2004T ): (2) feai et all tood \: (3) Ichung et all feoioh : (4) 



McDonald et all ( 1201 if ). 
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Table 2 





T TyiTi &v q n H T ,fwxT&v T ,1 in c 

KJ UUC1 dllLi JJUWCI U 11111 Lo 


r\T\ ( . c\ in tti T,n m it)/"ici'I~'1£*g 
wii v ' r i nil lift i el y ±j Liiiiiiivjoi ijico 




Cluster 


log(L .l-100Gev) 


log(Lo.l-lOOGev) 


log(Lo.l-lOOGev) 


Name 


(from L 8 _iooo Mm ; ±0.25 dex) 


(from L x 4 GHz ; ±0.19 dex) 






lower limit 


lower limit 


upper limit 


Coma (LF) a 


41.07 




42.55 


Coma 


40.50 


41.07 


42.55 


AWM7 


40.01 


40.17 


42.15 


Perseus 


39.73 


39.93 


43.55 


Hydra 


40.28 


40.55 


41.74 


A1367 


40.70 


41.19 


42.07 


Virgo 


40.04 


40.29 


41.61 


Bullet 


42.06 




44.66 


Ophiuchus 6 




40.66 


43.45 


A2029 6 




43.06 


43.46 


A2142 fe 




40.66 


43.54 



Notes. Luminosities are reported in erg s _1 . The uncertainties in the lower limit luminosities 
are from lAckermann et al.l (120121 ) and are a measure of the dispersion in the relationships 
used to calculate L .i_iooGeV- Upper limits on gamma-ray luminosity (column 4) assume a 
power law spectrum w ith a 7 = 2 in converting from the 0.2 — 100 GeV energy band, as in 
Ackermann et al.ll2010l . to 0.1 — 100 GeV for direct comparision with lower limits estimated 



from IR and radio luminosities. 

( a ) LF is luminosity function. See text for details. 

( b ) The data for these clusters is for BCGs only, not the full clusters. 



- 16 - 



45 



41 



Perseus 



'=c 43 
o 



12 



bC 
O 



41 



10 



39 



I 



AWM7 

\, Hydra 
Virgo T 



Coma Coma (LF) 

A1367 



II 



I 



I I 



I 



Bullet 



I 



10.5 



11.0 



11.5 12.0 
log(L8-1000jLtm) (i©) 



12.5 



13.0 



Fig. 1. — Gamma-ray luminosity versus total infrared luminosity. Blue points correspond 
to lower limits on gamm a-ray luminosity from star formation , calculated using r elatio ns in 
Ackermann et al.l (120121 ) . The uncertainties are reported in lAckermann et al.l ( 120121 ) and 
are a measure of the dispersio n in the Ln,i-inoOftV-Ls- -inorvrri relationship. Green arrows 



correspond to upper limits from lAckermann et al.l (120101 ) . 
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Fig. 2. — Gamma-ray luminosity versus radio luminosity. Blue points correspond to 
lower limits on gamma -ray luminosity from star formati on, calculated using rela tions in 



Ackermann et all ( 120121 ). The uncertainties are reported in lAckermann et all (120121 ) and are 
a measure of the dis persion in the Ln.i- i nnORV ~Li.4CtH7, relationship. Green arrows correspond 
to upper limits from lAckermann et all ( 1201 01 ). 



